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depending on parameters of the system. We study the resulting phase diagram in detail and inves- 
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mation and using a full real-space generalization of DMFT. 
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I. INTRODUCTION 



At the present moment, ultracold atoms in optical 
lattices can be considered as one of the most powerful 
tools for testing Hubbard- type models [1]. In these sys- 
tems parameters like the hopping amplitude, interaction 
strength, number and type of atomic species, lattice ge- 
ometry and dimensionality, can be tuned in a wide range, 
which makes it experimentally possible to verify predic- 
tions of theoretical models in different regimes. Despite 
incredible progress in the field of ultracold atoms, and, 
in particular, observation of short-range magnetic corre- 
lations in a recent experiment [2| , probing quantum mag- 
netic phenomena still remains a challenging goal for ex- 
perimentalists in this field, because the range of entropies 
(and temperatures) that are currently accessible is still 
several times higher than the upper limit for observing 
magnetic long-range order [3] . Once this important prob- 
lem has been solved, it is believed that ultracold atomic 
systems can give significant insight into high-temperature 
superconductivity [4[ and will be highly promising for 
quantum simulations [E[ in general. 

In this paper, we focus on a generalization of the Fermi- 
Hubbard model by considering two-component mixtures 
with repulsive interactions, which have an imbalance 
both in the hopping amplitude (corresponding to differ- 
ent effective masses) as well as in the population of the 
two components. In ultracold atomic mixtures, this gen- 
eralization was theoretically studied for the case of at- 
tractive interactions in the context of competing super- 
fluid and density- wave ground-states [6], Fulde-Ferrell- 
Larkin-Ovchinnikov superfluidity and other long-range 
ordered phases [7|, [8[ . For the case of repulsive interac- 
tions, it was investigated in the context of itinerant fer- 
romagnetism (Stoner instability) 0, [loj- Antiferromag- 
netic phases in ultracold imbalanced Fermi-Fermi mix- 
tures with moderately strong interactions were studied 
only separately in the case of population pj]-[l3| or mass 
imbalance HHi^. Evidently, so far the magnetic order- 



ing phenomena in the case, where both types of imbal- 
ance are present, have not been considered. Our paper 
aims at bridging this gap. 

The paper is organized as follows. In Sec. [XT] we in- 
troduce the model and derive an effective Hamiltonian, 
which allows us to interpret our numerical results. In 
Sec. IIIII we briefly outline the main steps in the general- 
ization of our numerical approach, dynamical mean-field 
theory (DMFT) [17[ , that are important for a proper ac- 
count of magnetic ordering effects in the case of mass and 
population imbalance. Sec. [TV] is devoted to our results 
and their discussion. 



II. MODEL 

We consider a Fermi-Hubbard Hamiltonian of the fol- 
lowing type: 

h =- J2 + hx -) + u Yl 

(ij) cr i 

+ ^2^2(Vi-^)7iia, (1) 

i <j 

where t a is the hopping amplitude of fermionic species 
a — {t?i}> cla (^icr) is the corresponding creation (anni- 
hilation) operator of species a at the lattice site z, the no- 
tation (ij) indicates a summation over nearest-neighbor 
sites, and U is the magnitude of the on-site repulsive 
(U > 0) interaction of the two different species with cor- 
responding densities and (n^ = c\ a Ci a ). In the 
last term, Vi is the external (e.g., harmonic) potential at 
lattice site z, and /i CT is the chemical potential of species 
a. Note that we have taken the harmonic potential to be 
independent of the atomic species. The Hamiltonian ([1]) 
implies a single-band approximation; in other words, we 
consider the case of a sufficiently strong lattice potential, 

Mat > 5£ r . 



It should be noted that within this model, the so-called 
mass imbalance depends only on the hopping amplitudes 
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t a , which for sufficiently deep lattices are given by [l8[ 

t a w ^^ ra ^/ 4 exp(-2v^v), (2) 
v 7r 

where ^ = ^f/^, ^ra = ^ 2 ^ 2 /2 ra a is the recoil en- 
ergy, k is the wave number determined by the wavelength 
of the laser forming the optical lattice, and m a is the mass 
of species a. The amplitude of the lattice potential 

can be different for the two components ^ ^iit^ 

which results in the possibility to realize an imbalance in 
the hopping amplitude even for different hyperfine states 
of the same atom (m^ = m^) [l9[. As for the population 
imbalance, its magnitude depends not only on the chemi- 
cal potentials /v, but also on other parameters, including 
the mass-imbalance At = (t^ — t±)/(t^ + t±). In order to 
characterize the population imbalance quantitatively, we 
introduce the polarization P = (N t - N i )/(N t + TV;), 
where N a is the total number of particles of type a in 
the system. 

It is important to mention that the introduced mass- 
imbalance parameter At can be experimentally tuned in 
a wide range. For example, according to Ref. [7], for a 
6 Li- 40 K mixture it can be effectively tuned from 0.3 to 
0.85 by varying the intensity (in the range of 1 Watt) and 
detuning (in the range of 2 nm from the magic wave- 
length) of the laser beams forming the optical lattice. 
Other systems where hopping imbalance can be realized 
and tuned in different ranges include the 171 Yb- 173 Yb 
mixture [20[ , spin-dependent optical lattices for homonu- 
clear mixtures and mixtures of alkali and alkaline-earth 
fermionic atoms. By approaching the limit of large im- 
balance, At — >• 1, these systems even allow for studies 
of the Falicov-Kimball model [21] that is used to model 
certain solid-state materials. 

By using the Schrieffer- Wolff transformation in the 
limit t(j <C U near half- filling, + « 1, we can map 
the Hamiltonian (pQ) to an effective spin model [22|, |23| . 
For the system under study, this transformation results 
in the anisotropic Heisenberg (XXZ) model, 

^eff = J\\ ^2 sfsf + J ± ^(sfs? + sYsJ) 

-Am£^' (3) 

i 

with coupling constants 

J,, =2(t 2 t + tl)/U, J±=^tJU, (4) 

and a chemical potential difference A/i = ji^ — ji^ that 
plays the role of an external magnetic field. Sf' (R = 
{X, Y, Z}) are usual spin- 1/2 operators on the lattice site 
z, Sf" = ^cla&apCiPi where a R are the Pauli matrices. 
It should be noted that here and below, all "magnetic" 
characteristics refer to a pseudospin made of two differ- 
ent species, thus the anisotropy in the spin model is not 
related to a spatial direction or quantization axis in the 
original lattice. 



The anisotropy in the Hamiltonian (j3]) originates from 
the imbalance in hopping amplitudes. This reduces the 
initial SU(2) rotational spin symmetry of the balanced 
mixture to the lower Z2 x U(l) spin symmetry, where 
the symmetry group Z2 corresponds to reflections about 
the XY plane [Ising type; first term in Eq. ([3])] and U(l) 
corresponds to spin rotations in the XY plane [second 
term in Eq. J3J)]. The symmetry of the model is further 
reduced to U(l) if one has a nonzero chemical potential 
difference A/i. 

According to Eqs. (j4]), we note that the coupling Jy is 
always larger than or equal to J±. Hence, from Eq. ([3]) 
one concludes that at A/i = and At ^ the ground 
state of the system is an easy-axis antiferromagnet (Ising- 
type; Z-AF). As was pointed out in Ref. [15], the ex- 
citation spectrum is gapped in this case and the Z-AF 
phase is also permitted in low dimensions (d < 3; as 
the Mermin- Wagner theorem applies only to continuous 
symmetries). Evidently, this also holds for a nonzero 
but small A/i < (Jy — J_l), when the system's ground 
state corresponds to a ferrimagnet (Z-AF with an addi- 
tional net magnetization in Z direction). In the opposite 
case, A/i 7^ and At = 0, the ground state is a canted 
antiferromagnet (easy-plane antiferromagnet with a net 
magnetization in Z direction), which obeys the Mermin- 
Wagner theorem and was studied in Ref. [13|- Therefore, 
in the region of intermediate imbalances, A/i ~ (J\\—J±), 
one should expect a phase transition between these two 
different types of magnetic ordering. 

The effective Hamiltonian ((3]) gives a good understand- 
ing of the types of ordered phases arising in the system 
under study. However, in order to have a more complete 
picture of the structure and quantitative characteristics 
of the magnetically-ordered phases arising in optical lat- 
tices at nonzero temperature and governed by the Hamil- 
tonian (jTJ) , it is necessary to use nonperturbative numeri- 
cal approaches. In this article, we apply dynamical mean- 
field theory, which is well suited for the description of 
long-range ordered phases in high-dimensional systems 
and is able to fully capture effects of inhomogeneity and 
finite-size that are usually present in optical lattices. 



III. METHOD 

Dynamical mean-field theory [I?} is an approach that 
allows to bridge two limits on the lattice: the nearly- free 
fermion gas and the strongly-interacting (atomic) limit. 
This method maps the lattice problem (which is, in gen- 
eral, intractable) to an impurity problem (which is chosen 
to be numerically solvable), thus substituting the full ac- 
tion by an effective one. Despite the fact that it is a non- 
perturbative approach, it is still an approximate method, 
since it treats the lattice self-energy as a local (which for 
homogeneous systems means: momentum-independent) 
quantity, thus neglecting nonlocal quantum fluctuations. 
DMFT becomes exact in the limit of infinite dimensions, 
d = 00 (i.e., large coordination number z ^> 1). Although 
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it is not an exact method in the case of square and cubic 
lattice geometries (z = 4 and z = 6, respectively), re- 
sults obtained by DMFT can be used as a reference point 
both for experiments and for more sophisticated meth- 
ods, such as quantum Monte Carlo simulations, which 
could be computationally rather demanding due to the 
presence of a sign problem in the case under study (for 
recent results, possibilities and limitations see [24] and 
references therein). 



A. Impurity model and solver 



Most commonly, for solving the auxiliary impurity 
problem in DMFT, the lattice model (pQ) is mapped onto 
an Anderson impurity model (AIM). The model Hamil- 
tonian contains the full local physics of the lattice prob- 
lem, but the non-local terms are represented by non- 
interacting bath degree's of freedom. In the case under 
study this Hamiltonian can be written in the following 
form: 



IM 



EE 6 ^^ + ^^Vi <T (a] (T d <T + h.c.) 

1=2 a 1=2 cr 

n s n s 

E E A 4>- + E E w iMJ-° + h - c -) 

1 = 2 cr 1 = 2 cr 



where a and its opposite a represent the spin indices, 
a = {t,l}, the index I = {2, . . . , n s } labels the number 
of the bath's orbital in AIM with n s being the cut-off 
number. In our calculations we use the exact diagonal- 
ization (ED) solver [25| with n s = 5, such that there are 
4 orbitals. The two terms in the first line of Eq. ([5]) cor- 
respond to the energies of electrons in the bath and the 
hybridization between the bath and the impurity, respec- 
tively. The second line represents the anomalous terms 
that are important for obtaining quantities correspond- 
ing to magnetic ordering in the XY plane. The terms in 
the last line have a direct correspondence to the Hamilto- 
nian (pp), except the last (auxiliary) term, which is used 
in calculations as a small initial perturbation to break 
the remaining U(l) symmetry of the spin model (j3]). 

Using the standard technique (see, e.g., Ref. [26|) of 
eliminating bath degrees of freedom in the effective action 
corresponding to the Hamiltonian (|5]), we find analytical 
expressions for the Weiss Green's functions, which rep- 
resent effective dynamical fields acting on the impurity 



site: 

n s 

Q~ X {iUn) = iu n + /V - E [ V ll( iuJ n ~ 

1=2 

+2Vl ff W la A l +W? a (iLJ n -e l<r j\ , 

n s 

1=2 

HVitVn + Wi t Wi^)Ai + V n W n (iuj n - e n )} , (6) 

where K\ = (iuj n - e^)(iuj n - £4) - A^, uj n = 7r(2n + 
1)//? is the Matsubara frequency and j3 is the inverse 
temperature, /3 = 1/T (we use units such that fc# = 1). 

Within the ED solver, the basis states of the finite- 
dimensional Hilbert space are given by 



t t 

^1 j n 2i ' 



> 2 , . 



n s I 



(7) 



with rip = 0, 1 and ^2, p n° = n a . Note that the anoma- 
lous terms in Eq. ([5]) mix the sectors and (i.e., 
the magnetization s z is not conserved), which therefore 
cannot be diagonalized independently. Although the to- 



tal charge n 



is still conserved, this leads to 
a significant increase of the numerical effort in diagonal- 
ization and subsequent calculations of the corresponding 
Green's functions. At finite temperature, these are calcu- 
lated from the full set of eigenstates \i) (with eigenvalues 
Ei) according to 



G aia2 (iuj n ) =^E 
(e~ 



-fSEj 



Z 



(i\lAM\j)(3\dl\i) 



^O"! ^CT 1 

Ei — En 



(8) 



(9) 



where Z = Yli e /3Ei ls the partition function. Next, 
following Ref. |27j. the self-energies can be defined as 



Saa = U 



E T(jG<7<7 E acr G uu 

Gn - - _ r 1 ^ ' 
aa^crcr ^ ctct 

E ycfG <j<j E a(J G era 

Gq-Q-Gq-Q- G fjjj 



(10) 

(11) 



In practice, we solve the impurity problem by ob- 
taining the quantities (l5])-(fTT|) for given parameters U , 
fi a and f3 from the original lattice problem (pp) and 
for a particular set of auxiliary Anderson parameters 
{ziaiVia^Wia, A^}, which is updated in each DMFT it- 
eration. The self-energies (JTJJJ) and (fTTj) then allow us 
to calculate the Green's functions corresponding to the 
initial lattice problem. Below we consider two main ap- 
proaches for evaluating these Green's functions: i) two- 
sublattice DMFT, which is important for obtaining phase 
diagrams for homogeneous (infinite) systems with mag- 
netic order and can also be used in combination with a 
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local- density approximation (LDA) to analyse trapped 
gases, and ii) real-space DMFT, which describes finite 
inhomogeneneous (trapped) systems without further ap- 
proximations. 



B. Two-sublattice DMFT 

Bipartite structures, such as states with antiferromag- 
netic order, can be described in an appropriate way by 
introducing two sublattices. Within the DMFT approach 
one then needs to solve the impurity problem twice on 
two adjacent sites of the original lattice. It is impor- 
tant to note that here, in contrast to the case of bal- 
anced mixtures, the observables corresponding to differ- 
ent sublattices (denoted below by s = 1,2) cannot be 
directly associated with observables corresponding to dif- 
ferent species (denoted by a -l)- Hence, we define the 
lattice Green's functions in the following (generalized) 
way: 



CT1CT2 



(iu> n ) = f deD{e)[A-\e)l% 

J — Z 



(12) 



with 



/ c w 



and 



A*) 
As) 

Scrcr 



tt 
/-(I) 

Hi 

—t^e 
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Hf 

AV 
Hi 





—t^e 
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A 2 ) 
Hi 



o 

A*) 
Ht 

A 2 ) 

Hi 



(13) 



(14) 

(15) 



where z denotes the lattice coordination number {z = 4 
and z = 6 for square and cubic lattices, respectively), 
D(e) is the normalized density of states, f* deD(e) = 1, 
the explicit form of which is known for a particular lattice 
geometry. The local self-energies appearing in Eqs. ([T4]) 
and JTSJ) are taken from the impurity solver [see Eqs. ([TO]) 
and ([IT])]. 

To complete the self-consistency equations of the 
DMFT scheme, we define the Weiss Green's functions 
from the Dyson equation, 



1 

criCT2 



[G (,) W^ + SlN. (16) 



where G^ s ^ is a 2 x 2 block of the matrix ([T2]h the in- 
verse of which is taken separately for different sublattices 
s = 1,2. By using the obtained Weiss Green's functions 
in the minimization procedure [according to Eqs. (j6j) and 
applying a conjugate gradient method], we find a new set 
of Anderson parameters, which is then used in a subse- 
quent DMFT iteration. These iterations are performed 
until final convergence, i.e., until the initial and final 
Weiss Green's functions coincide. 



When an inhomogeneous system is considered, the 
two-sublattice DMFT introduced above can be used 
in combination with LDA. The main advantage of 
LDA+DMFT is that this approach allows to consider 
large systems in 3d, as the numerical effort scales approx- 
imately linearly with the system size for axial-symmetric 
trapping potentials. The drawback of this approach that 
it fails to reproduce the detailed structure close to the 
boundaries of the ordered phases, i.e., it does not ac- 
count for a possible proximity effect. The mentioned ef- 
fect can be accounted for within another generalization: 
the real-space dynamical mean-field theory, which was 
firstly introduced in Refs. [1|[29|. 



C. Real-space DMFT 

The main idea of real-space DMFT (R-DMFT) is not 
to divide the lattice problem into several sublattices, but 
to solve the impurity problem on each lattice site corre- 
sponding to the original finite-size system directly. Then, 
after we obtain the self-energies E c v^ a - 2 (za; n ) [see Eqs. ([TO]) 
and ([TT]) ] for each lattice site i = {l,...,iV}, we collect 
them in the real-space matrix consisting of inverse local 
Green's functions and hopping elements, 



G 



/ Ctt 

^(i) 
Hi 

h 





V . 



\ 



AD 
Hi 



h 
o 

o 



A*) 

nt 

A 2 ) 

Hi 

h 
o 



a(2) 

Hi 


h 



a(3) 

u t 

/-(3) 

Hi 



.(3) 
ii 



(17) 



/ 



which is hermitian and of size 2A^ x 2A^. Here, the diag- 
onal matrix elements C^l have the form: 

(&(i"n) = iu n + H*-Vi- £$(«<;„), (18) 

and the off-diagonal elements (i l )(ioj n ) are defined ac- 
cordingly to Eq. ([T5]) . 

As in the two-sublattice case, we close the DMFT-loop 
with the lattice Dyson equation, 

[G^Hicvn)]-^ = [G«(io, n )]-V a + ^Jl 2 (iio n ), (19) 

where is a 2 x 2 block of the real-space Green's 

function matrix G obtained by inversion of ([TT]) . Finally, 
one defines a new set of Anderson parameters (as in the 
two-sublattice case) for each lattice site. 

In case of a large system size (when the total number of 
lattice sites N > 10 3 ) the inversion of the real-space ma- 
trix ([T7]) becomes a time-consuming task in the numerical 
calculations (in comparison with the impurity solver with 
a moderate number of bath orbitals). Nevertheless, even 
with this limited total number of lattice sites R-DMFT 
is capable of a proper description of proximity-induced 
effects in lattice systems with magnetic order. 
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IV. RESULTS 
A. Homogeneous systems 

1. Unpolarized mixtures with mass imbalance 

Firstly, let us discuss the effects originating only from 
the mass imbalance in the system. We set fi^ = /i^, 
which results in balanced populations (P = 0) in a ho- 
mogeneous system. According to Sec. [Ill we note that 
the magnetic ground state of this system at half filling is 
a Z-antiferromagnet (Z-AF) for any nonzero value of the 
mass imbalance. 

A first important effect that we want to emphasize is 
that, according to mean-field analysis and existing Monte 
Carlo calculations in the limiting cases At = and At = 
1 (see Refs. |30U3lj|. respectively), the critical (Neel) tem- 
perature increases with the mass imbalance. Indeed, by 
taking the mean-field definition for the Neel tempera- 
ture in the Heisenberg model [32], T/y = 6 JS(S + l)/3, 
where S is the fermions' spin, one obtains for constant 
U/t, where t = (t t + ^)/2, and J = Jn [see Eq. (0}] 



T N (At) = (l + At 2 )T N {0). 



(20) 



In order to prove that this effect also takes place for the 
Hubbard model (pQ) at moderately strong interactions, 
U ^ we numerically calculated T c by means of two- 
sublattice DMFT for a cubic lattice. The results are 
shown in Fig. DJa), which clearly confirms this behav- 
ior. According to Ref. [l5[, this phenomenon can be ex- 
plained as caused by the suppression of quantum fluctu- 
ations (due to the emergence of an energy gap for one 
of the Goldstone modes) in systems with nonzero mass 
imbalance. 

Another interesting effect concerns the entropy analy- 
sis for homogeneous systems. As pointed out in Ref. pjl, 
one can approach much closer the ordered state for mass- 
imbalanced mixtures with the same entropy values, com- 
pared to the balanced system. In order to show this, we 
present in Fig. DJb) the dependence of the minimal tem- 
perature that can be reached for constant entropy when 
minimizing over the interaction strength U/J (thus as- 
suming an adiabatic change of the interaction strength) 
as a function of the mass- imbalance parameter. Note that 
for s = 0.7 this curve has a maximum. This is because for 
small mass-imbalance the lowest temperature is reached 
for large U/t such that the system is in a Mott insulator. 
In this regime the minimal temperature increases when 
the mass-imbalance is made larger. For even larger mass- 
imbalance the lowest temperature is reached at small U/t 
such that the system is in the Fermi liquid phase. In this 
phase the lowest temperature decreases with the mass- 
imbalance, giving rise to the observed maximum. 

It is important to note that it is now well-established 
that DMFT quantitatively overestimates the critical tem- 
perature and critical entropy values in comparison with 
more exact methods for balanced mixtures. According 




0.4 0.6 

Mass Imbalance At 

FIG. 1. (a) Critical temperature for Neel ordering in the half 
filled Hubbard model (/it 4 — U/2) versus mass imbalance at 
different values of the interaction strength in a cubic lattice, 
obtained within DMFT. (b) Dependence of the minimal tem- 
perature that can be reached for a fixed entropy s per particle 
by minimizing over the interaction strength U/t, as a function 
of the mass- imbalance parameter At. 



to our analysis, this is also the case for systems with 
nonzero mass imbalance. However, we should stress that 
the DMFT approach becomes more precise in the limit 
of large mass imbalance not only regarding the quanti- 
tative estimates of the critical temperature, but also in 
the critical entropy, the value of which does not depend 
on mass imbalance within the dynamical mean-field de- 
scription for U > zt, s™ f = In 2 « 0.69. This can be 
seen, in particular, from the comparison with known re- 
sults for the Heisenberg model [33[ (U ^> t, At —> 0), 
where s c « 0.34, and results for the Ising model [34j 
(U ^> t, At —> 1), where s c ~ 0.56. Hence, the discussed 
advantages of imbalanced mixtures should be even more 
pronounced in studies based on more exact methods. 

It should be noted that in addition to the spin-density 
wave in mass-imbalanced mixtures one also observes a 
weak charge-density wave (CDW) modulation. This is 
directly related to the presence of both Neel ordering 
and mass imbalance in the system and corresponds to 
the fact that in the Z-AF state the sites occupied by a 
heavier component have an enhanced double-occupancy 
due to hopping from adjacent sites, which are occupied 
by the lighter component (in the same way, the oppo- 
site mechanism works for the sites occupied by a lighter 
component). The magnitude of this CDW, according to 
the estimates presented in Ref. [35| for U > zt, is pro- 
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FIG. 2. Phase diagram for a population- and mass- 
imbalanced mixture at half filling and finite temperature, 
obtained within DMFT for a cubic lattice. Dashed and 
solid lines correspond to the first- and second-order phase 
transition lines, respectively. The CDW parameter is de- 
fined in terms of the double-occupancy Di = (n^n^), c = 
(A - A+i)/(A + A+i); miy = (S^ x ) - (S?g). 



portional to At • (U/t)~ 2 . Hence, the CDW is more pro- 
nounced at moderate interaction strength, and it van- 
ishes in the large U/t limit. 



2. Polarized mixtures with mass imbalance 

The presence of both population and mass imbalance 
in ultracold mixtures according to the effective Hamilto- 
nian (j3]) results in competition between different types 
of antiferromagnetic ordering. In Fig. [2] we present 
the phase diagram showing the structure of the ordered 
phases in the intermediate parameter region at half fill- 
ing [for moderate polarizations P < 0.8 and U > zt, the 
condition of half filling is fulfilled by taking the average 
chemical potential ft = U/2, where ft = (/i^ + /i^)/2, 
but for larger population imbalances it must be adjusted 
by hand] . When the hopping amplitude of both species 
is low in comparison with the interaction strength, the 
system is in a paramagnetic Mott insulating state. In 
the opposite case, the system is in the unordered Fermi 
liquid phase, which is compressible and characterized by 
an enhanced double-occupancy of the lattice sites. In 
the central region, two different phases appear: a canted 
antiferromagnet, which is characterized by a staggered 
magnetization in X direction and a net magnetization in 
Z direction, and a ferrimagnet (labelled by Z-AF in this 
diagram), which has both a staggered and a net mag- 
netization in Z direction. Note that in the canted-AF 
phase, the occupation of each site by heavy and light 
species is equal, thus an additional CDW doesn't ap- 
pear in this case. In the region of intermediate mass 
imbalance, a first-order phase transition (with a narrow 
coexistence region) takes place between the phases with 




0.4 0.8 1.2 6t t /f7 0.4 0.8 1.2 Qt t /U 



FIG. 3. Set of phase diagrams for a population- and mass- 
imbalanced mixture at half filling, T = 0.4t and different A/i, 
obtained within DMFT for a cubic lattice, including easy- 
axis antiferromagnetic (Z-AF), charge-density wave (CDW), 
ferrimagnet ic, canted AF, unordered Mott insulator (PM-I) 
and unordered Fermi liquid (PM-FL) phases. 



different magnetic order. 

We now discuss how the phase diagram presented in 
Fig. [2] changes when A/i (which defines the polariza- 
tion P) is increased (see also Fig. [3]). At A/i = the 
whole region with magnetic order corresponds to the Z- 
AF phase, except the central part (namely, the diagonal 
line = ti), where the AF order can be set in any di- 
rection. When population imbalance appears, the cen- 
tral phase characterized by canted order develops. This 
phase increases in width (towards the regions with large 
mass imbalance) when A/i becomes larger. As for the 
ferrimagnetic region, it shrinks not only from the side 
of the canted-AF phase, but also from the sides of un- 
ordered phases, as the net magnetization in Z direction 
suppresses weak Z-AF order in the region close to phase 
boundaries. At large population imbalances (A/i > 1.3t 
for T = OAt) the ferrimagnetic phase completely van- 
ishes, leaving only the region with canted-AF order in 
the central part. 

The antiferromagnetic order is therefore at a given 
temperature most unstable against imbalances in the 
triple-point region of a diagram of the type presented 
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FIG. 4. Dependence of critical polarization destroying AF 
order on the mass- imbalance parameter, obtained by DMFT. 



in Figs. [2] and [3j where canted, ferrimagnetic and para- 
magnetic phases coexist. In Fig. |4] we show this effect 
for different temperatures of the imbalanced mixtures. 
The critical polarization clearly reveals the sensitivity 
of the magnetically-ordered phases to an increase of im- 
balance in the triple-point regions, where corresponding 
minima are observed. Also, the effect that ferrimagnetic 
states are suppressed by population imbalance can be 
seen from this plot by comparing the left-hand (canted) 
and right-hand (ferrimagnetic) sides of curves presented 
in Fig. |U Interestingly, DMFT predicts that at zero tem- 
perature antiferromagnetic order develops for any polar- 
ization. This is different for the case of attractive inter- 
actions and superfluidity, which is destroyed at a finite 
Chandrasekhar-Clogston limit for the polarization [36l — 
[38[, which was beautifully demonstrated in experiments 
with ultracold Fermi gases [39|. 



B. Finite systems in a harmonic trap 

Finally, we turn our analysis to the case when an ad- 
ditional external trapping potential is present in the sys- 
tem. For simplicity, below we assume that the trap- 
ping potential Vi in Eq. (pQ) has an axial-symmetric form, 
Vi = Vorf/a 2 , where Vq is the strength of the harmonic 
potential, vi is the distance from the lattice site i to the 
trap center, and a is the lattice constant. Evidently, all 
the results of this section can be extended to include 
anisotropics usually present in real experiments. 

We start our analysis by applying the local-density ap- 
proximation in combination with two-sublattice DMFT 
introduced in Sec. IIIIBI Within LDA, we perform cal- 
culations for a particular lattice site by taking y$ = 
Mcr°^ + Vi in the DMFT calculations. In all results pre- 
sented in this section the average chemical potential in 
the center of the trap is taken as (/i!^ + /ij°^)/2 = [7/2, 
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FIG. 5. Real-space density and magnetization distribu- 
tions for a cubic lattice in a harmonic trap obtained by 
LDA+DMFT. Parameters used in (a)-(d): U = 12*, T = 0.2*, 
A* = 0.2, and V = 0.1*(n/a) 2 , where a is the lattice con- 
stant. 



such that the system is at half filling in the trap center 
[at least for \A/jl\ <C U (A/i = /i!^ — which is the 

most relevant case for our investigation. 

Performing calculations for fixed mass imbalance and 
varying chemical potential differences A/i in the trap cen- 
ter, we obtain distributions of the total filling and mag- 
netization in various directions, which are presented in 
Fig. [5l There are several details in those plots that are 
worth discussing. First of all, we note that for nonzero 
mass imbalance and A/i = one has a (globally) polar- 
ized mixture due to a "ferromagnetic" shell originating 
from a wider distribution of the lighter component in 
a trap. Hence, for the purpose of obtaining Z-AF or- 
der in the trap center, one should account for this effect 
and adjust the total polarization (or make the harmonic 
trap also species dependent). One could also use this 
adjustment in a combination with the additional cool- 
ing mechanism in mass-imbalanced mixtures, pointed out 
in Ref. [l6[ (removing from the system a part of the 
heavy component that carries larger amount of entropy 
at U = 0). Note however, that as long as one only speci- 
fies the chemical potential difference, there is no problem 
here. 

Possible advantages of the Z-AF phase could be not 
only the issues related to cooling as pointed out in 
Ref. [15] and in Sec. IIV A 11 but also a more convenient 
detection, e.g., by single-site resolution imaging with a 
quantum gas microscope [40[ , where the alternating den- 
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sity of one spin (species) component can be directly de- 
tected. In contrast, in the canted-AF phase, the density 
of the spin components does not alternate from site to 
site, thus one has to analyze the behavior of the double- 
occupancy [4lL l42j|. merge nearest-neighbor lattice sites 
[H or introduce additional methods to detect nearest- 
neighbor spin-spin correlation functions in the XY plane, 
such as Bragg spectroscopy (43j or noise correlations anal- 
ysis nn . 

Is is worth noting that the charge-density wave struc- 
ture peculiar to the Z-AF phase is clearly seen in the dis- 
tributions of the total particle number in Figs. [5fb),(c) 
(double lines in the central part), while it is absent in 
Figs. Efa),(d), where the canted-AF phase developed 
in the bulk. As one sees, the shell structure also has 
an interesting dependence on the population imbalance. 
The magnetization in it increases with total polarization 
(adding light and removing heavy particles from the sys- 
tem), but with the decrease of total polarization the "fer- 
romagnetic" shell related to mass imbalance does not 
vanish. Instead, according to Fig. [5fa), a double-shell 
structure develops with inner and outer "ferromagnetic" 
shells originating from population and mass imbalances, 
respectively. 

There are also discontinuities present in the magneti- 
zation close to the boundaries of the antiferromagnetic 
phases: since the applied LDA+DMFT approach does 
not include proximity effects these sharp features are not 
smeared out by the trap. Moreover, in these regions this 
approach has a rather bad convergence. In Fig. [6] we 
show that by accounting for the proximity effect within 
R-DMFT, the magnetization in trapped imbalanced mix- 
tures has a smooth behavior with a larger region of sta- 
bility of the antiferromagnetic phase than predicted by 
LDA. 



V. CONCLUSIONS 

We studied antiferromagnetically-ordered phases that 
emerge in two-component ultracold fermionic mixtures 
with mass- and population imbalance. Our analysis was 
based on dynamical mean-field theory and its real-space 
generalization at finite temperature. 

It is pointed out that two types of imbalance favor dif- 
ferent types of antiferromagnetically-ordered phases: the 
ferrimagnetic phase is favored by mass imbalance, while 
the canted-antiferromagnetic phase is favored by popula- 
tion imbalance. In the absence of population imbalance, 
we demonstrated within DMFT the advantages of mass- 
imbalanced mixtures, i.e., an increase of the critical tem- 
perature and the possibility to cool the system by adia- 
batic change of the interaction strength to lower tempera- 
tures than possible for balanced mixtures. Guided by the 
exact values of the critical entropy in the limiting cases of 
the Heisenberg and Ising model, we argued that DMFT 
gives a better prediction for the critical entropy for imbal- 
anced mixtures compared to balanced ones, which makes 
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systems with mass imbalance even more advantageous for 
the purpose of observing magnetic ordering phenomena 
in optical lattices. 

In the presence of both mass and population imbal- 
ance we obtained the finite-temperature phase diagram 
with the corresponding first-order phase transition be- 
tween the different AF states. We revealed that AF or- 
der is most unstable against thermal fluctuations in the 
triple-point regions. To this end, we performed a stability 
analysis of the ordered phases against population imbal- 
ance at different temperatures. At zero temperature we 
found that for all polarizations antiferromagnetic order 
develops. 

We also obtained real-space density and magnetization 
distributions of imbalanced mixtures in a harmonic trap. 
It is shown that, depending on the total polarization, 
the mass-imbalanced mixture can have different ordered 
phases in the bulk and different magnetic shell structures. 
The detailed description of these effects could help not 
only in preparing the mixture closely to its equilibrium 
state, but also in the detection of antiferromagnetic cor- 
relations in ultracold gases. 
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